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Abstract. We investigate the static critical behaviour of a uniaxial magnetic layer, 
with finite thickness L in one direction, yet infinitely extended in the remaining 
d dimensions. The magnetic dipole-dipole interaction is taken into account. We 
apply a variant of Wilson's momentum shell renormalisation group approach to 
describe the crossover between the critical behaviour of the 3-D Ising, 2-d Ising, 3-D 
uniaxial dipolar, and the 2-d uniaxial dipolar universality classes. The corresponding 
renormalisation group fixed points are in addition to different effective dimensionalities 
characterised by distinct analytic structures of the propagator, and are consequently 
associated with varying upper critical dimensions. While the limiting cases can be 
discussed by means of dimensional e expansions with respect to the appropriate upper 
critical dimensions, respectively, the crossover features must be addressed in terms of 
the renormalisation group flow trajectories at fixed dimensionality d. 



PACS numbers: 05.70.Jk, 05.10.Cc, 75.40.Cx 
1. Introduction 

Layered magnetic systems are of obvious technological importance, e.g., in magnetic 
storage devices. Presently, magneto-optic storage devices based on amorphous GdTbFe 
and TbFeCo alloys are in commercial use, allowing for storage densities of up to 10 8 
bits/cm 2 . In the future, however, Pt/Co, Pd/Co, and Co/Au multilayers of a few 
Angstrom thickness are expected to lead to four times as high storage capacities. On 
the purely theoretical side, they furthermore represent an intriguing model system 
with the competing exchange- and dipolar-dominated universality classes on the one 
hand, and the dimensional crossover between three and two effective dimensions on 
the other hand. The mathematical description of the associated crossover features 
within the framework of the renormalisation group theory is technically challenging, as it 
involves both varying analytic structures of the propagators, and different upper critical 
dimensions for the associated renormalisation group fixed points. The aim of the present 
paper is to demonstrate that such a unifying crossover description is feasible, even 
with the aforementioned difficulties at hand. For a layered uniaxial magnetic system 
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of finite width L in one direction and with dipolar interactions taken into account, the 
renormalisation program is explicitly performed to one-loop order for the entire crossover 
regime by means of a variant of Wilson's momentum shell method at fixed transverse 
dimension d. Moreover, it is shown that the scheme described here encompasses all the 
relevant universality classes involved, and the critical exponents at each renormalisation 
group fixed point may be computed within an e expansion about the respective upper 
critical dimensions. 

In the theory of phase transitions there has been long-standing interest in the 
critical properties of finite systems as well as in the influence of the magnetic dipole- 
dipole interaction on static and dynamic critical behaviour. The first investigations for 
the crossover behaviour between the different critical regimes in dipolar systems date 
back thirty years: In their pioneering work, Larkin and Khmel'nitskii [|IJ showed that the 
upper critical dimension of a uniaxial magnet with magnetic dipole-dipole interaction 
is reduced by one as compared to the pure Ising system. The mathematical description 
of the crossover in a three-dimensional system with easy axis and dipolar forces 
therefore necessitated a renormalisation scheme capable of dealing with varying upper 
critical dimensions associated with the different non-trivial fixed points. Modifying 
the generalised dimensional regularisation method by Amit and Goldschmidt 0, such 
a description was developed in Ref. || for the crossover between the Ising and 3- 
d uniaxial dipolar fixed points of a uniaxial dipolar magnet, and in Ref. [|J for the 
crossover from isotropic to directed percolation. This approach was then applied to an 
even more complex scenario by Ried et al. [||], who considered the crossovers between 
all four possible fixed points of a uniaxial dipolar system, taking the finiteness of the 
anisotropy into account. It should be noted, though, that all these calculations employed 
the harmonic dipolar propagator originally derived in Ref. Yet this propagator is 
adequate only in the case of a three-dimensional infinite system, or, more generally, 
for any system of the same dimensionality as the one used for computing the effective 
magnetic dipole-dipole interaction. It is therefore not applicable for 3-d electrodynamics, 
if the system itself is two-dimensional, for example, or is finite in one or more directions. 

Static critical behaviour in finite systems with purely short-range interactions was 
investigated in Refs. [7j ||, ||, Here, the asymptotic critical behaviour was addressed, 
and no crossover scaling functions were calculated. Dimensional crossovers were studied 
in the context of quantum phase transitions by Chakravarty, Halperin, and Nelson for 
quantum antiferromagnets [O], and by O'Connor and Stephens for an Ising system 



between four and three dimensions ||12|| . Yet we are not aware of any work concerning 
the critical behaviour in a finite system that takes into account the magnetic dipole- 
dipole interaction. 

In this paper we go beyond the cited literature in describing the crossover scenario 
for a uniaxial magnetic system with dipolar interactions that is infinitely extended in 
d dimensions, but has a finite thickness L in one additional direction. For simplicity, 
we assume that the anisotropy energy is larger than all other relevant energy scales, 
and hence only take into account the Ising-type up-down symmetry. In the case 
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of d = 2, such a system is supposed to exhibit four different regimes of non-trivial 
critical behaviour, in addition to the classical one described by mean-field critical 
exponents; these are described by the 3-D Ising, 2-d Ising, 3-D uniaxial dipolar, and 2-d 
uniaxial dipolar universality classes, respectively. As will be shown, these regimes are 
characterised in addition to varying effective system dimensions by different analytic 
forms of the propagator, which originate in the dipole-dipole interaction. Consequently, 
three different upper critical dimensions d c emerge: For the Ising systems, D c = 4 and 
d c = 4, respectively, for the 3-D uniaxial dipolar system D c = 3, and finally the upper 
critical dimension for the truly asymptotic critical behaviour of a 2-d uniaxial dipolar 
system with an easy axis in the system's plane will be demonstrated to be d c — ~. 

This paper is organised as follows: In section ||] we introduce the effective 
hamiltonian (free energy functional) and specifically the propagator (two-point 
correlation function in Gaussian approximation) that is to be used for our calculations. 
In section [| we explain our variant of the momentum shell renormalisation method 
devised to mathematically describe the crossover scenarios. In section f| the resulting 
renormalisation group flow equations are derived to one-loop order and discussed. 
Section |5| contains the resulting limiting cases within the e expansions in the vicinities 
of the appropriate upper critical dimensions. In section ^ the crossovers between the 
critical fixed points are discussed in terms of the renormalisation group trajectories. We 
conclude with a summary and brief discussion in section [?|. 



2. Effective hamiltonian 



At first we need to deduce an effective free energy functional that is capable of capturing 
the crossover features we are aiming at. However, at this point we encounter the 
difficulty that the analytic forms of the propagators (Gaussian two-point correlation 
functions) describing uniaxial dipolar magnetic systems are quite different depending 
on the dimensionality D = d+1 of the system. We employ periodic boundary conditions 
to avoid surface effects, in which we are not specifically interested here. Rather we are 
concerned with the effect of the system's limited extent in one space direction as such. 

Denoting the part of the wave vector k in the direction of the easy axis of the 
magnetisation as q, and the perpendicular component as p, k = q + p, the propagators 
look as follows: In the case of a three-dimensional system, the relevant part of the inverse 
propagator is of the form Go(q,p) _1 = r + ck 2 + gq 2 /p 2 ||13|| , while for a two-dimensional 
system containing the easy axis, Go(9 ) p)~ 1 — r + ck 2 + gq 2 /k, with k = \k\ [[14]]. On the 
other hand, if the easy axis is aligned perpendicular to the 2-d plane, the result is of the 



form Go(q,p) = r + ck — g k \14\. Here r = T — T c denotes the temperature distance 



to the mean-field transition temperature T c °, c is a constant describing the stiffness of the 
system with respect to inhomogeneous order parameter configurations, and g measures 
the strength of the dipole-dipole interaction as compared to the exchange interaction. 

As our task is the description of a system with finite thickness, we have to apply 
some care in the derivation of the appropriate form of the propagator. The contribution 
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to the free energy stemming from the dipole-dipole interaction between the magnetic 
moments connected to the spin variables Sr at lattice sites R reads 

W D = g £ Y,^(Ri,R2)S« Rl S^, (1) 

Ri^R 2 a/3 

with the dipole tensor in real space being defined as 

A {RuR2) ~ IR.-R^ \R^M 5 ' ( ) 

Here, the indices a and (3 denote real space vector components. In Fourier space, the 
dipole sum for a cubic lattice becomes 

A^(fc) = Af d { P ) + — e%? E (p ~ f )a(p 7 9)P e ~ Hlp ~ g| > ( 3 ) 

a o p ^o g \p-g\ 

if the spins are oriented in the plane, and 

A^(k) = A^(p)--J2 eiqP J2\P-9\e- Wp - 9] , (4) 

if the spins are oriented perpendicular to the plane |14[]. A%H(p) represents the dipole 
tensor for the two-dimensional system, with wave vectors p and associated 2-d reciprocal 
lattice vectors g. p is the discrete coordinate of a lattice plane along the finite direction, 
and ao and q indicate the corresponding lattice constant and wave numbers, respectively. 
For the later treatment by means of the renormalisation group technique, we need only 
retain the relevant long-wavelength contributions, and may thus simplify the above 
expressions by restricting the reciprocal wave vector sums to g = only. This part 
yields the correct limits for the propagators of both the three- and the two-dimensional 
system. 

Finally, we expand the resulting expressions to obtain the long-wavelength limits. 
For the easy axis lying in the plane we find the following form: 

a 3 ad 

A^(q,p) = (1 + h) **- + f 2 — q 2 + / 3 pV 

p p 

+ Up a p P P + hp a p P q 2 + 0(p 3 q 2 ,p\pY) . (5) 

Denoting with 2N + 1 the number of lattice planes in the finite direction, L = 2N ao, 
the iV-dependent functions fa read: 

fx = 2N , (6) 

f2= -^N(l + 3N + 2N 2 ) , (7) 

h = -a N(l + N) , (8) 

f 4 = ^N(l + 3N + 2N 2 ) , (9) 

h = ^N 2 (l + N) 2 . (10) 
If the easy axis is aligned perpendicular to the plane, we arrive at 

A^(q,p) = -(1 + A)p- f 3 p 2 - f 2 q 2 p + 0(p 3 q 2 } p 4 ) (11) 
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instead, with the iV-dependent functions 
/i = 2N , 

„2 

f 2 = -^N(l + 3N + 2N 2 
6 v 

h= - ao N(l + N) . 



(12) 
(13) 
(14) 



We note that the long-wavelength expansion does not commute with taking the limit 
of infinite system thickness, N — > oo. Consequently, the above formulas are strictly 
valid only for finite system thickness L. Nonetheless we shall see that the method to be 
described below will yield the correct results even for the three-dimensional system. 

We are now in the position to write down the effective hamiltonians that serve as 
starting points for our renormalisation group treatment. Dealing with a finite system 
in one direction, we symbolically write the sum over all Fourier modes as 



(27T) 



/p I £ ... , 

0<|p|<A Li 



n e Z . 



(15) 



q =^n 



In the limit of the fully infinite system, L — > oo, this definition is to be replaced with 

/ . . . = — / d d p I dq , where D = d + 1 (16) 
Jk [2ti) u jo<|p|<a J 

denotes the dimension of the entire system. In both situations, A represents a cutoff 

whose numerical value is not relevant for the critical behaviour. Therefore we shall 

use A = 1 throughout the paper (in appropriate dimensionless units). Notice that we 

only apply the cutoff for the wave vector components in the plane, and not for the 

perpendicular component. 

Thus the final effective free energy for a layer with easy axis in the plane assumes 

the form 



n 



2d.\ 



1 

2Jk 



+ 



ro + Cok 2 + g{1 + fl) d + gf2 d q 2 

p p 

+ ghp 2 i +gfap'iP + gf5Plq 2 ] s (k)s (-k) 

5 (fci)5 (fc 2 )So(fc 3 )5 (fc 4 ) *fx>) 



(17) 

4! Jk\ Jk-2 Jk?, J k4 

where q again is the wave vector component perpendicular to the system, and p the 
in-plane part, with p\ denoting the component pointing along the direction of the easy 
axis. The hamiltonian for a layer with easy axis perpendicular to the plane reads 

1 



H 



2d,. 



r + c k - g(l + fi)p- gf 3 p - gfoq P\ S (k)S (-k) 



4! Jk\ Jkz JkiJkn 

For a three-dimensional system we find instead [1 



\i=l 



n 



3c/ 



1 

2 Jk 



r + c k 2 + g 



P 2 



S (k)S (-k) 



+ 



Uo 

4! J ki Jk2 Jk% Jk4 



5o(fci)5o(fc 2 )So(fc 3 )5 , o(fc4; 



ii k: ) 



(19) 
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where q denotes the component of k along the direction of the easy axis, and p is the 
part of k perpendicular to the easy axis. 

Thus, henceforth we are concerned with the description of critical phenomena with 
those three different hamiltonians. The free energy of the three-dimensional system 
cannot be recovered any more from the layered hamiltonian (fTBD by taking the limit of 
infinite layer thickness iV — > oo. Nevertheless we will show that all the limiting cases, 
including the three-dimensional uniaxial dipolar one, are described correctly within our 
scheme. 

3. Renormalisation scheme 

As already mentioned, the expected critical regimes are clearly characterised by different 
system dimensions, as well as different analytical structures of the propagators, leading 
to different upper critical dimensions. Hence we need a renormalisation scheme capable 
of dealing with these different fixed points in a unified manner, and sufficiently flexible 
as to describe the crossovers between them. In addition, we aim to capture the 
crossover features with the correct transition temperature T c (L,g) (see appendix A) of 
the finite system, and not with that of the infinite system as in Refs. f?|, [| || [[(J , where 
furthermore the studied systems were infinite in at most one direction (and consequently 
do not display a proper phase transition). As our system asymptotically becomes two- 
dimensional, we have to circumvent problems with infrared divergences appearing in 
field-theoretic approaches, e.g., dimensional regularisation, for such low-dimensional 
systems [0, [| [12], [OJ. One more aspect is that we must include the full gradient terms 
in order to describe the crossover between Ising and uniaxial dipolar behaviour, and 
cannot restrict ourselves to merely the most relevant parts in one of the limits as in 
Ref. |B|. 

In order to accomodate all these requirements, we apply a variant of Wilson's 
momentum shell renormalisation method, similar to the one used previously by 
Chakravarty, Halperin, and Nelson for the crossover theory of a quantum non-linear 
sigma model [fLl|j . (Certainly to one-loop order, the ensuing wave vector shell integrals 
over p, followed by the summations over q, are technically considerably easier to 
perform in this approach as compared to the field theory version.) One renormalisation 
transformation now consists of two steps: Starting from the effective hamiltonian 
functional 7i, we first integrate out the short-wavelength order parameter fluctuations 
in the partition function Z(7i) in the d directions of the plane, and then sum over the 
wave vector components perpendicular to the plane (without any cutoff). In wave vector 
space, the modes to be integrated out are depicted in figure [3]. In the second step, the 
wave vectors k, spin variables S k , and system thickness L are rescaled with an as yet 
arbitrary scale factor b > 1, and a spin rescaling factor ((b) which eventually has to be 
chosen appropriately in terms of b: 

k! =bk , (20) 

S' k , = ((b)- 1 S k , (21) 
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Figure 1. The modes being integrated out during one renormalisation step. 



U =b' 1 L . (22) 

After this transformation, we arrive at a new free energy functional Ti! of the same form 
as before, but with renormalised coefficients fulfilling Z(7i) = Z(H'). By successively 
evaluating partial partition functions, and subsequently applying scale transformations, 
we approach the interesting infrared behaviour of our theory, and in effect map 
the perturbationally inaccessible critical regime onto a region in parameter space 
where perturbation theory may be applied. As usual the critical regions themselves, 
characterised by algebraic singularities and scale-invariance, are represented by fixed 
points of the renormalisation (semi-)group transformation. 

We apply this scheme in the framework of the one-loop approximation, however 
refraining from any e expansion, for in our problem there is no unique and common 
upper critical dimension. We choose our rescaling factors isotropic in contrast to, e.g., 
Ref. ||16|| ; this induces apparently diverging parameters under the renormalisation flow 



HI. Yet these may be readily absorbed into appropriately defined effective couplings, 
by means of which we can describe the associated crossover features 0, H, H, |5|- 
contrast to Ref. we define the spin rescaling factor ((b) such that the gradient term 
of the hamiltonian functional remains unchanged during the renormalisation procedure. 
Again, this helps us capturing the entire crossover regime, while the advantages of the 
scheme in Ref. [[I7j], where the term with the lowest wave vector power is held fixed, are 
gathered through the construction and discussion of suitable effective couplings. 
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4. Renormalisation group flow equations 

In this section, we compute the differential renormalisation group flow equations for a 
magnetic layer with the easy axis in the plane, a layer with easy axis perpendicular to 
the plane, and a three-dimensional system. Our starting points will be the different 
hamiltonians fll7D , ([Tp, and (|19|), respectively, for these systems. Following through 



the renormalisation scheme as described in the previous section above, with b = e , 
5 —>■ 0, for the finite systems we arrive at the following differential flow equations for 
the system thickness L, the coefficient c of the gradient term, and the strength g of the 
dipole-dipole interaction: 

%=-L, (23, 

Before moving on to the flow equations specific to the different systems, we discuss 
these differential equations common to both finite systems. The flow equation 
for the layer thickness L is exact, and is solved by L = Lqq~ 1 with t — > +oo in the 
infrared limit. Thus, the fixed points for the layer thickness are L = oo and L = 0, 
with L = being infrared-stable. Hence the asymptotic critical behaviour is that of 
the two-dimensional system: As the correlation length diverges, the system thickness 
becomes irrelevant in the renormalisation group sense (as indicated by the negative sign 
on the r.h.s. in equation (|23|) .) The coefficient c of the gradient term in the free energy 
functional is held fixed during the renormalisation procedure through the choice of the 
spin rescaling factor (, whence equation (|24"D follows from definition. The strength of 
the non-analytic dipole-dipole interaction g cannot be renormalised perturbationally, 
and consequently equation (^) simply follows from dimensional analysis and holds to 
all orders in perturbation theory. Its solution reads g = goe e , and the fixed points for 
the dipole strength are g = and g = oo, and the latter is infrared-stable. As to be 
expected, the critical behaviour of the system is eventually dominated by the dipole 
interaction. This is a consequence of the long-range character of the dipolar forces in 
contrast to the short-ranged exchange interactions. For the derivation of the remaining 
flow equations we have to distinguish between the different systems under consideration. 

First we focus on the layer with the easy axis in the plane. The remaining differential 
flow equations read 

f = 0, (26) 



d£ 
d£ 
d£ 



-2/ 2 , (27) 
- h , (28) 
-2/ 4 , (29) 
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Figure 2. Definitions of pi and the angle 9. 



dfa _ 

IF' " 3/5 



following, to one-loop order, simply from dimensional analysis, and 



(30) 



dr 
di 

du 
di 



u 

2r + - 



A 2 r 27r 

— / de 



2 (2vr) d L Jo 



L/2 



coth I — . 



a(8)b(6) 



a(0) 
2\b{&) 



(31) 



A 2 



2tt 



L/4 



sinh(§ 



6(0)- 



X 



cosh I — A 



q(0) 



a(0) 



2\ &(#) oi , , l , / A /ni^i 



(32) 



to one-loop order. Here indicates the angle between the wave vector component in the 
plane p and the direction of the easy axis p\ (see figure ||). The ^-dependent functions 
a and b are defined as follows: 

a(0) = r + cA 2 + g(l + fa) A cos 2 6 + ^/ 3 A 2 cos 2 6 + ^/ 4 A 3 cos 2 6 , (33) 
6(0) = c + #/ 2 A cos 2 + ^M 2 cos 2 . (34) 

While the coefficient fi remains invariant under the renormalisation group flow, the 
irrelevant parameters fa, fa, fa trivially flow towards their stable fixed points fa = fa = 
/4 = 0. Notice that there exist no fixed points with fa = oo, because the functions fa 
were only defined for a finite system thickness. The more intricate flow equations for the 
temperature variable r and the non-linear coupling u will be discussed in the following 
sections. 

For a layer with easy axis perpendicular to the plane, on the other hand, we find 
the following flow equations 

dfa 

de 

dfa 
di 







-2/ 2 



(35) 
(36) 
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di 

from power counting, and 



^="/s. (37) 



Here the functions a and 6 are defined as follows 

a = r + (c - g~h) A 2 - g(l + J\) A , (40) 
b = c-gf 2 A-gf 4 A 2 , (41) 

and we have introduced the used the dimension-dependent factor 

A d 

K = 7 — r-: — t . (42 

Closer inspection reveals that in the case of an easy axis perpendicular to the layer and 
with infinite anisotropy, the propagator is not positive definite. Inevitably therefore, 
the terms in the sums in equations ( |3~8D and (|39"D will become negative under the 
renormalisation group flow. This indicates an instability at non-zero wave vector, 
and the ground state of the system becomes inhomogeneous, with the characteristic 
wavelength of the emerging spatial structures given by that instability wave vector (the 
typical size of the structures appearing in the low-temperature phase is estimated in 
appendix B). At any rate, it is impossible to analyse the asymptotic critical behaviour 
with the effective hamiltonian flT8|) which was obtained via an expansion about a 
presumed homogeneous ground state. Hence, in the remainder of this paper we shall 
concentrate on the study of the critical behaviour of a layer with in-plane easy axis, and 
on the three-dimensional limit. 

For a three-dimensional system, finally we obtain the renormalisation group flow 
equations 

£=<>, (43) 



d^ 
dg 



2<7, (44) 



— = 2r + - (45) 



4 



(c+&) (r + cA 2 ) 



^- = {A-D)u-\u 2 K . (46) 

di 8 ^(c+^) (r + cA 2 ) 3 

These equations do not appear as the limit L — > oo of the flow equations of the layer 
(pip and (|52"|), where the effective free energy has already been adjusted to the specific 
geometries. In order to arrive at the correct 3-D expression, one has to start from the 
correct, full hamiltonian (|19|) , as the long- wavelength limit and the limit — > oo do 
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not commute, as already mentioned in section 0. Again, c is held fixed by construction, 
and the renormalisation group flow for the relevant dipole-dipole interaction has the 
unstable fixed point g = and the infrared-stable fixed point g = oo. 

Before discussing the different critical regions and the crossover between them 
within our scheme, we focus on the obvious limiting cases of the general flow 
equations, and compute the critical exponents within dimensional e expansions about 
the appropriate upper critical dimensions, respectively. 



5. Limiting cases in e expansions 

5.1. Three-dimensional I sing limit 

We arrive at the limit of the flow equations at the unstable fixed point g = and L = oo 
by expanding the hyperbolic functions in the differential flow equations ([31]) and ([32]) 
for the layer in the limit of large argument. Assuringly, we can alternatively take the 
flow equations for the three-dimensional system (|45|) and (f46|) with g = 0. With either 
approach, we obtain 

dr u K . . 

-T7 =2r + - 47 
^ 4 v / c ( r + cA 2) 

du 3 2 K 

— — — €U u 2 . , (48) 

d£ 8 J c ( r + C A2)3 

where e = 4 — D gives the distance to the upper critical dimension D c = 4 of the system. 
From these equations we get the critical exponents at the non-trivial fixed point to first 
order in e, 

r* = - e -ck 2 + 0{t 2 ) , (49) 

Ku ' = hfk? +0{ * 2) - (50) 

as usual via the eigenvalues of the linearised flow equations. The results, consistently 
expanded to first order in e = 4 — D are 
1 / e 



2v 1 + 6 + ° (e2) J ' (51) 
rj = + 0(e 2 ) , (52) 

1= l + i + 0(e 2 ), (53) 
o 

in accordance with the well-known results for the single- component 4 model (see, e.g., 
Ref. |15|), with e = linD = d + l = 3 dimensions. 



5.2. Three-dimensional uniaxial dipolar limit 



Next, we discuss the second limiting case in three dimensions, namely the uniaxial 
dipolar fixed point at g = oo and L = oo. We have to start from equations ([45]) and 
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(pj|). Upon introducing the effective expansion parameter || 

v e s = u/y/g , (54) 
we arrive at the flow equations 

d^eff 3 2 K A 

-M= £V «-S V °" (r + cA^/2 ' (56) 
where now e = 3 — -D was used for the deviation from the upper critical dimension 
D c = 3. In .D = d + 1 = 3 dimensions, we thus expect the mean-field critical exponents 
with logarithmic corrections. This is in agreement with previous work addressing the 
critical behaviour of a three-dimensional uniaxial dipolar magnet |], |3|. ||. 

5. 3. Two-dimensional Ising limit 

In the limit of absent dipole-dipole interaction g = and zero thickness L = we have 



to expand the hyperbolic functions in the flow equations (|3l|) and ( j32|) for vanishing 
arguments. With the effective coupling 

v eS = u/L, (57) 

this procedure gives the flow equations 

*L - 2r I VeR - (58) 

d^eflf 3 2 K 

-M="«-2 V *ff (r + cA^r (59) 

here, however, e = 4 — c2 (= 2inZ) = c2+l = 3 dimensions). Linearising in the vicinity 
of the non-trivial fixed point 

r* = - JcA 2 + 0(e 2 ) , (60) 
o 

Kv: s = 2 -ec 2 A* + 0(e 2 ) } (61) 
we find the critical exponents to order e = 4 — d(= 2 in two dimensions): 

v = \ (l + ~ + ^(^)) , (62) 

rj = + C(e 2 ) (63) 

7 = l + J + 0(e 2 ). (64) 
6 

Thus, this low-dimensional limit is also correctly included in the general flow equations. 
Note that in terms of the expansion parameter e, these exponents of course coincide with 
the three-dimensional Ising case above. Remember, however, that there e = 4 — D = 1 
and here e = 4 — d = 2in the relevant physical situation D — d + 1 = 3. 
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5-4- Two-dimensional uniaxial dipolar limit, in-plane easy axis 

This limit is especially interesting for two reasons: First, it represents the asymptotic 
critical behaviour of the system under consideration; and second it is intriguing on its 
own right because we are aware of only one other reference which addresses the critical 
behaviour of a two-dimensional uniaxial dipolar system with the easy axis in the plane 



Yet, in Ref. |LJ| the critical exponents are not presented explicitly. Using the 
same approximation for the hyperbolic functions as before and introducing the effective 
expansion parameter 

v eS = — ^— , (65) 

we obtain the asymptotic flow equations 

o> „ _Wcff A 2 

M r + 2<W ^( r + cA 2 ) (1 + /i) A ' 1 J 

dVeff ~ f7 (67) 



d£ \2 J n 2 2 d 



7T 



r + cA 2 ) 3 (l + /!)A 



Consequently, we identify the upper critical dimension here as d c = 7/2. The non-trivial 
fixed point to order e = \ — d is 

r * = _i cA 2 + C(e 2 ) , (68) 
3 

< ff = [(1 + A)A] 1/2 (cA 2 ) 3/2 + 0(e>) . (69) 

Linearisation near this fixed point leads to the following asymptotic critical exponents 
to first order in e = | — d 

" = \ (l + I + 0(e^ , (70) 
V = + O(e 2 ), (71) 
l = l+ e - + 0{e 2 ) . (72) 

Once again, these exponents look superficially identical to those of the short-range Ising 
model; yet the critical dimension and the value of e is different: in d = 2, e = 1/2. In 
summary, all the asymptotic limits are correctly included in the general flow equations. 



6. Crossover analysis 



In the previous section, we have demonstrated the consistency of our general approach 
with the e expansions near the respective critical dimensions of the limiting universality 
classes. We now proceed to the application of our scheme to the description of the actual 
crossovers between the different critical regions. 
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6.1. General procedure 

At first we want to describe the general procedure to be used. We use the calculated 
renormalisation group flow equations to investigate the scale dependence of the effective 
interaction parameters. Specifically, we map the critical region onto a region in 
parameter space that is perturbationally accessible. In the infrared region, we 
can, e.g., evaluate the static susceptibility perturbationally, and thus obtain the 
static susceptibility x m t ne critical region via the solution of the renormalisation 
group equation that takes the form (the Fisher exponent rj vanishes in the one-loop 
approximation) 

X (r ,g ,u ,L ,k) = e 2£ X (r(i),g(i), u(£), L(£), k(£)) . (73) 

Therefrom we can infer the effective exponent 

. , dhiY -1 . , 

= in^r < 74 » 



as function of reduced temperature r = T ■ This method was used first in |T9| in the 



framework of an e expansion; but here we make no use of any dimensional expansion, 
but perform all all our calculations numerically at a fixed dimension. The one-loop 
recursion relations are then solved without recourse to any additional approximation. 

It is vital here to introduce an appropriate effective expansion parameter that yields 
the correct 3-D and 2-d limits as given in the previous section. This effective coupling has 
to have a finite limit at the asymptotic fixed point below the upper critical dimension 
of the system. As suggested by the flow equations in section f|, we use the effective 
couplings 



/■*r 1 L 

v e ff = u I dO = coth — A 



(75) 



a(0) 

'° ^Ja(9)b(9) \ b(6) 

for the layer with in-plane easy axis, and 

2?r u 
Vr + cA 2 y/g 

for the three-dimensional case. Then within our methodical framework, the 
critical regions are again represented by fixed points of the renormalisation group 
transformation, and we may obtain the critical exponents from the eigenvalues of the 
linearised transformation in the vicinity of the fixed points (yet without applying any 
further approximation like dimensional expansions). 



6.2. Limiting cases at fixed dimension 

The appearance of renormalisation group fixed points to first order in an expansion does 
not necessarily guarantee the existence of fixed points if we refrain from an e expansion. 
At first, therefore, we have to re-calculate the non-trivial fixed points and the associated 
critical exponents at fixed dimensionality D = d + 1. The actual fixed points thus 
obtained are listed in table |l|. The eigenvalues of the linearised flow equations can be 
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Table 1. Rcnormalisation group fixed points at fixed dimension to one-loop order. 
Fixed point r* v* s 

3-D Ising -|(4 

2- d Ising -(4 - 

3- D uniaxial dipolar 
2-d uniaxial dipolar — (~ 



Table 2. Numerically determined fixed point values and critical exponents. 



Fixed point 


r* 




L* 


7 


V 


3-D Ising 


-0.25 


78.95 


OO 


1 


1/2 


2-d Ising 


-0.25 


78.95 





1 


1/2 


3-D uniaxial dipolar 








00 


1 


1/2 


2-d uniaxial dipolar 


-0.33 


105.27 





0.94 


0.47 



- D)cA° 



2 + |(4 



D) 



(4-D) 



cA 2 



- d)cA 2 /(4 - d + 6) 

- d) cA 2 / (l-d + 3) ff (f - d) cA 2 / (|-d + 3 



16tt 

±g(4-d)cA 2 /(4-d + 6) 




easily calculated analytically, and are not presented explicitly here. In table || we give 
the numerical values of the fixed points as well as the critical exponents v and 7 for the 
system with two infinite dimensions (D = d + 1 = 3). As is readily seen, the numerical 
results for the critical exponents cannot be judged as good approximations to their 
realistic values. The reason is very likely the inadequacy of the one-loop approximation. 
An unfortunate artefact of the system under consideration is that the actual fixed point 
values as well as the critical exponents for the 3-D Ising and the 2-d Ising limit turn out 
identical to this order. This happens to be the case only in this dimension D = 3, as 
can be seen from figure [3], where we plot the numerical values of the critical exponent 
v at both the higher and lower dimensional fixed points, respectively, vs the dimension. 
This artificial coincidence only appears at the dimensions d = | and d = 2, the physical 
dimension in which we are interested. Nevertheless we shall argue for the usefulness of 
our scheme in describing different crossover features. 

6.3. Dimensional crossover 

The first crossover that we want to analyse within our scheme is the dimensional 
crossover for an Ising layer (dipolar strength g — 0). In order to circumvent the 
unfortunate exponent coincidence mentioned above, we evaluate the crossover features 
between 3.2 and 2.2 dimensions. In this situation, we are not only able to calculate 
the flow equations describing the crossover, but also the effective critical exponent 
7eff as a function of the reduced temperature r. The corresponding renormalisation 
group fixed point values are listed in table [| The effective exponent 7 e fr describing 
the crossover is calculated numerically for different layer thicknesses L, and the results 
are plotted in figure £|. It can be seen that at a temperature sufficiently far from the 
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0.52 



""2.0 3.0 4.0 

d 

Figure 3. Ising fixed point values of the critical exponent v as function of dimension. 
Table 3. Numerical fixed point values for an Ising layer with d = 2.2. 



fixed point 


r* 


Ks 


L* 


V 


3.2-D Ising 


-0.2105 


96.0267 


00 


0.5090 


2.2-d Ising 


-0.2308 


105.2601 





0.5070 



transition temperature, the susceptibility exponent takes on the value of the higher- 
dimensional fixed point, then with decreasing temperature goes through a minimum 
which is obviously independent of the layer thickness, and finally reaches the asymptotic 
value of the lower- dimensional fixed point. The thicker the layer, the system remains 
the longer in the higher-dimensional critical region as the temperature is lowered, and 
reaches the later the asymptotic critical region. This is fully to be expected based on 
straightforward considerations. 

6-4- Crossover exchange- dominated / dipole- dominated 

Next we apply our scheme to the fully three-dimensional system, where the crossover 
occurs between the region with dominating exchange interactions and the asymptotic 
region with dominating dipole-dipole interaction. Because the numerical values of the 
critical exponents happen to be identical at both these fixed points within our scheme 
at fixed dimension D = 3, see table the crossover cannot be noticed at all in the 
temperature dependence of the effective critical exponent j e ff. However, we can discuss 
the crossover by means of the flow diagrams in the r — v e g plane. In figure [|, the flow 
diagram is depicted for different initial temperatures (which cannot be distinguished in 
this graph, however). The arrow indicates the starting point of the trajectories. These 
first approach the unstable Ising fixed point (filled circle on the right), before finally 
approaching the asymptotically stable fixed point that represents the uniaxial dipolar 



d-dim. FP 
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10" 



Figure 4. Layer with vanishing dipole-dipole interaction: Effective critical exponent 
7 e ff vs reduced temperature r ~ r — r c for systems of different thicknesses: N = 20 
(full line), N = 50 (dashed), TV = 100 (dotted), TV = 500 (long-dashed). 



regime (filled circle on the left). Eventually the trajectories diverge in the temperature 
variable r due to the numerically inevitably finite distance from the critical point (even 
small inaccuracies in determining the true critical temperature will eventually blow up 
under renormalisation). In figure ^| we depict the same scenario for a stronger dipole- 
dipole interaction. In this situation, the unstable fixed point is not approached as closely 
as before, and the stable fixed point is approached more closely before the trajectories 
eventually diverge to r — > oo due to the finite distance from the critical temperature. 

6. 5. Crossover features of the dipolar layer 

At last, we address the description of the crossover scenario for a dipolar magnetic layer 
with in-plane easy axis. As the starting point is the hamiltonian fllTD, which is only 
defined for finite system thickness, we cannot reach the fixed point corresponding to the 
three-dimensional uniaxial dipolar region. The crossover will thus occur between the 
critical regions of the three- and two-dimensional Ising systems, and the two-dimensional 
uniaxial dipolar regime, which represents the asymptotic universality class. As the 
numerical fixed point values of the effective coupling, the temperature variable, and 
the effective critical exponent of the static susceptibility 7 e g- of the two- and the three- 
dimensional Ising fixed point are identical, we illustrate the flow trajectories in three- 
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Figure 5. Infinite three-dimensional dipolar system: Flow of the effective coupling 
v e g and temperature variable r for different initial temperature values, with g = 1CP 5 
and Mo = 5: r w r c (full line), r — r c = 6.42 ■ 10~ n (dotted), r — r c = 6.66 ■ 10~ 9 
(dashed), r — r c = 1.67 • (long-dashed). 



dimensional parameter space, consisting of r, t> e ft, and the thickness-dependent variable 
-j^-j- as third coordinate. In terms of this quantity, the three- and two-dimensional Ising 
fixed points assume the value 1 and 0, respectively, and thus become distinguishable. 
In figures |7|, || and |9| we show the trajectories with initial values in the vicinity of 
the critical temperature for increasing values of the dipole strength. At first the 
trajectories approach the three-dimensional Ising fixed point, then the two-dimensional 
Ising fixed point, and finally the asymptotic two-dimensional uniaxial dipolar fixed point. 
Eventually even the tiny numerical deviation from the critical temperature gives rise to 
a divergence in r again. For stronger dipole interaction, the Ising fixed points are not 
approached as closely In figure [LCj we furthermore depict the temperature dependence 
of the trajectories: The higher the temperature the sooner the trajectories start to 
diverge. 
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Figure 6. Infinite three-dimensional dipolar system: Flow of the effective coupling 
v e g and temperature variable r for different initial temperature values, with go = 1CP 3 
and u = 5: r « r c (full line), r - r c = 4.44 • lCT 10 (dotted), r - r c = 4.04 • 1CT 8 
(dashed), r — r c = 7.40 • 10~ 7 (long-dashed). 

7. Summary 

In this paper we have analysed crossover scenarios for a uniaxial dipolar magnetic layer. 
We have derived the appropriate effective free energy functional in the long- wavelength 
limit for a finite system with thickness L in one spatial direction; the resulting dipolar 
hamiltonian turned out to be considerably different from the one for the infinite three- 
dimensional system. In the long- wavelength limit, which does not commute with the 
limit L — > oo, the propagator of the three-dimensional system cannot be recovered from 
the propagator of the layer in the limit of infinite layer thickness. 

We have presented a unified renormalisation scheme, within which the critical 
behaviour of the layer with finite thickness can nevertheless be treated as well as 
that of the infinite system. This renormalisation procedure is required to be capable 
of capturing the crossovers between critical regions with different system dimensions, 
different analytical structures of the propagator due to the dipole-dipole interaction, 




Dimensional crossover in dipolar magnetic layers 




Figure 7. Dipolar layer with easy in-plane axis: Flow diagram for g 
10, r ~ r C) iV = 100. 




Figure 8. Dipolar layer with easy in-plane axis: Flow diagram for g 
10, r - r c ,iV = 100. 
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Figure 9. Dipolar layer with easy in-plane axis: Flow diagram for g = 10 5 ,u n = 
10, r ~ r C) iV = 100. 




Figure 10. Dipolar layer with easy in-plane axis: Flow diagram for g — 10 6 ,u = 
10, r -r c = 8.33-10- 10 ,2.33-10- 10 ,0.53-10~ 10 ,r - r c , (from left to right), 7V = 100. 
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and consequently different upper critical dimensions. Hence we had to refrain from 
any dimensional e expansion near any of the fixed points, and have derived the 
renormalisation group flow equations in the one-loop approximation at fixed dimension. 
Furthermore, we have introduced a common effective expansion parameter, as suggested 
by the flow equations, which leads to finite non-trival fixed points char act arising the 
different critical regions. As a test of the consistency of our approach, the critical 
exponents in the limiting cases of the 3-d Ising, 2-d Ising, 3-d uniaxial dipolar, and 
2-d uniaxial (in-plane easy axis) dipolar universality classes were calculated within e 
expansions with respect to the appropriate upper critical dimensions, respectively, and 
reproduced correctly known results from the literature. The critical exponents for the 
asymptotic limit of a two-dimensional uniaxial dipolar system with in-plane easy axis 
were, to our knowledge, computed explicitly for the first time by means of an expansion 
near the upper critical dimension d c = 7/2. 

Performing the calculations without the e expansion, however, reveals an 
unfortunate insufficiency of the one-loop approximation: The numerical values of the 
temperature variable, the effective coupling, and the critical exponents at the three- 
dimensional and the two-dimensional Ising fixed points happen to be identical in D = 3 
dimensions. In order to test our scheme, we have first evaluated the dimensional 
crossover of the effective critical exponent 7 e g of an Ising system from D = 3.2 to 
D = 2.2 dimensions. Second, the crossover features of the three-dimensional uniaxial 
dipolar system as well as the uniaxial dipolar layer were analysed and discussed in 
form of renormalisation group flow diagrams. Thus, we hope to have demonstrated 
the feasibility of the description of such very complex crossover features within the 
framework of our renormalisation group approach, which should be further applicable 
to a diverse range of interesting physical systems. 
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Appendix A. Critical temperature 

A number of the above calculations depended on the knowledge of the phase transition 
temperature which, due to thermal fluctuations, is shifted downwards as compared to 
the mean-field transition temperature. The correct location of the critical point (to one- 
loop order) was determined numerically through identifying that starting point of the 
temperature variable r, which separates those renormalisation group flow trajectories 
that eventually flow to r = +oo from the ones that eventually flow to r = — oo. 



Figures [Al| and [A2| show the dependence of the numerically obtained phase 



transition temperature r c upon the strength of the dipole-dipole interaction g and the 
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Figure Al. Phase transition temperature r c (g) vs dipole strength (No = 100, uo 
0.1). 



layer thickness N. As one would expect physically, r c grows both with increasing dipole 
strength and increasing layer thickness. 



Appendix B. Layer with out-of-plane easy axis 

As already mentioned, the propagator for the dipolar layer with out-of-plane easy axis, 
G (q, P)- 1 = r + c k 2 - g(l + fi)p- gf 3 p 2 - gfc q 2 V , (B.l) 

is not positive definite under the renormalisation group flow. This indicates an 
instability of the homogeneous ground state at some finite wave vector. Therefore 
this propagator cannot be used in a straightforward manner for the calculation of the 
critical behaviour of the system. Nevertheless we can draw some information regarding 
the instability from it. Thus, from the equation 
d 



dp 



q=0 



(r + ck 2 - g(l + /i) p - gf 3 p 2 - gf 2 q 2 pj =0 



(B.2) 



we may infer the typical size of the structures in the low-temperature phase, as function 
of the dipole strength g and the layer thickness L. The result is 



-^domain 



2tt 


c + g± 




2a y 




9 




-) 

«0 J 





(B.3) 



I.e., starting from a constant domain size at vanishing layer thickness, the domain 
size decreases with increasing thickness, runs through a minimum, and finally grows 
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Figure A2. 

0.1). 



Phase transition temperature r c (N) vs layer thickness (go = 1-10 5 , uq 



again linearly with L. 



The result that a layer with out-of-plane easy axis should 

; but the increase of 



condense inhomogeneously is in accordance with Refs. |20], [21|, |_ 
the domain size is predicted to be proportional to the square root of the layer thickness 

L 1/3 in Ref. |l] . 



in Refs. 20, 221, and 
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